Kaon B parameter from quenched domain-wall QCD * 
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We report on a calculation of Bk with domain wall fermion action in quenched QCD. Simulations are made 
with a renormalization group improved gauge action at /3 = 2.6 and 2.9 corresponding to a -1 w 2GeV and 3GeV. 
Effects due to finite fifth dimensional size N5 and finite spatial size are examined in detail. Matching to the 
continuum operator is made perturbatively at one loop order. We obtain Bk{^ = 2GeV) = 0.5746(61), where 
the error is statistical only, as an estimate of the continuum value in the MS scheme with naive dimensional 
regularization. This value is smaller but consistent with Bk(^ = 2GeV) = 0.628(42) obtained by the JLQCD 
Collaboration using the Kogut-Susskind quark action. Results for strange quark mass are also reported. 



1. Introduction 

The Kaon B parameter Bk is an important 
quantity to pin down the CKM matrix, thereby 
gaining an understanding on CP violation in the 
Standard Model. 

A crucial ingredient in a precision calculation 
of Bk is chiral symmetry. For this reason, the 
best result so far has been obtained H with the 
Kogut-Susskind (KS) quark action, which keeps 
U(l) subgroup of chiral symmetry at finite lattice 
spacings. Carrying out a systematic and exten- 
sive set of simulations to control 0(a 2 ) scaling 
violation and Ota- — ) errors that arise with the 

MS 

use of one-loop perturbative renormalization fac- 
tors, the value B K (fj, = 2GeV) = 0.628(42) was 
obtained in the continuum limit in the MS scheme 
with naive dimensional regularization (NDR). 
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The Wilson fermion action has a problem of 
explicit chiral symmetry breaking, which causes 
a nontrivial mixing between four-quark operators 
with different chiralities. This problem has been 
treated with several non-perturbative renormal- 
ization methods The resultant values of 
Bk are consistent with that of the KS action, 
but the numerical error is large. 

The domain wall fermion formalism offers a 
possibility of a calculation preserving full chiral 
symmetry ||. Here we report a summary of a 
quenched calculation of Bk with the Shamir's 
variant |?]] of the formalism and a renormaliza- 
tion group (RG) improved gauge action 0. The 
latter choice is motivated by our result || that 
chiral symmetry is much better realized with this 
action than with the plaquette gauge action. We 
may also expect that scaling behavior of Bk is 
improved with the use of the RG-improved ac- 
tion. We also report on the strange quark mass 
obtained from meson mass measurements in our 
simulation. 
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Table 1 



Simulation parameters together with the number 
of configurations analyzed shown in bold letters. 





(3 = 2.6 




(3 = 2.9 


a- L (GeV) 


1.81(4) 




2.87(7) 


N t 


40 




60 




16 24 


32 


24 32 


N 5 = 16 


122 76 


25 


76 50 


N 5 = 32 


50 







2. Numerical simulation 

Simulations are made at two values of coupling, 
(3 = 2.6 and 2.9 corresponding to a lattice spacing 
a" 1 = 1.81(4) GeV and 2.87(7) GeV determined 
by p meson mass, in order to check the scaling 
behavior. 

The choice of lattice size N% x N t x N$ is made 
as follows: (i) Our main runs use 24 3 x 40 x 16 
at /3 = 2.6, and 32 3 x 60 x 16 at [3 = 2.9. (ii) 
Spatial size dependence is examined at (3 = 2.6 
varying the spatial size from N a = 24 to N a = 
16 and 32: They correspond to the physical size 
aN a ~ 1.7,2.6,3.4 fm. (iii) Dependence on the 
fifth dimensional length N$ is also examined at 
(3 = 2.6, N a = 24 using N 5 = 16 and N 5 = 32. 

We choose m f = 0.01,0.02,0.03,0.04 for the 
bare quark mass, covering the range mps /my= 
0.4 - 0.8. The u-d and s quarks are assumed 
degenerate in the analysis for Bk- The domain 
wall height is taken to be M = 1.8. The number 
of configurations analyzed is given in Table |. 

We measure the Kaon B parameter, 

B _ (K\s^(l- l5 )ds^(l- l5 )d\K) 
K ~ § (K \sj^ 5 d\ 0) (0 \s^ 5 d\ K) ' 

and the matrix element divided by the pseudo 
scalar density, 

o _ (K |5 7m (1- 75)^(1 -!5)d\K) 
P (K\s l5 d\0)(0\s To d\K) [) 

which should vanish at m, — » 0. The extrac- 
tion of these quantities from Kaon Green func- 
tions follows the standard procedure (see, e.g., 
Ref. [|]]) using the Dirichlet boundary condition 
in time and wall sources for quark propagators. 
We simultaneously evaluate hadron masses. The 



physical point for quark masses is calculated by 
linearly fitting the light hadron masses m 2 PS and 
my as a function oi mf, and using the experi- 
mental values of m T /m p , mx/mp as input. 

3. Operator matching 

We carry out matching of the lattice and con- 
tinuum operators at a scale q* = l/a using one- 
loop perturbation theory |l(| and the MS scheme 
with NDR in the continuum. In the domain 
wall formalism the renormalization factor of an 
n-quark operator O n has a generic form 

Z = (1 - w 2 Q ) n / 2 ZZ /2 Z 0n (3) 

where wo = 1 — M and Z UKOn = 1 + 0(g 2 ). We 
apply tadpole improvement by explicitly moving 
the one-loop correction to the domain wall height 
M from Z w to wq , and by factoring out a tadpole 
factor u n l 2 — P n / S with P the plaquette from 
Zo n ■ A mean-field estimate appropriate for the 
RG-improved action is used for calculating the 
coupling constant 5 , |^g(/ i ) pTp - The continuum 
value at a physical scale e.g., fj, — 2.0 GeV, is ob- 
tained via a renormalization group running from 
q* = l/a to /a. 

For Bk the factor (1 — Wq) 2 Z^ cancels out, and 
the one-loop value given by the ratio 

Z Qi = l-(41og g *a+13.6)g 2 /(16^) 

Z\ (l-6.26 ff 2 /(167r 2 )) 2 1 > 

turned out to be very near unity i.e., Z^ s (q* = 
l/a) = 0.984(/3 = 2.6) and 0.988(/3 = 2.9). The 
Z factor at the scale p = 2GeV obtained with 
a 2-loop running becomes Zjf s (fi — 2GeV) = 
0.979(^ = 2.6) and 1.006(/3 = 2.9). 

4. Results for Bk 

4.1. Chiral property 

We start by investigating the chiral property of 
the system through Bp. The renormalized value 
of Bp, extrapolated to mj —Q linearly in mf, is 
plotted as a function of the inverse of spatial size 
in Fig. 0. 

Since results for the fifth dimensional size -/V5 = 
16 (filled circle) and 32 (open square) agree with 
each other, the observed non-zero value, albeit 
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[3=2.6, N 5 =16 
□ P=2.6, N s =32 
-» 0=2.9, N s =16 
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Figure 1. The matrix element Bp at rrif = as 
a function of spatial size. 
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Figure 3. Renormalized Bk as a function of spa- 
tial size. 



• B K at physical quark mass 



Figure 2. Bare Bk interpolated as a function of 
m f a at j3 = 2.6 for a 24 3 x 40 x 16 lattice. 

very small, should be due to a finite spatial-size 
effect. The decrease of the magnitude of Bp to- 
ward larger spatial sizes supports this interpreta- 
tion. The existence of finite-size effect is consis- 
tent with the result for the pion mass in Ref . |j] , 
where a non-zero pion mass of a magnitude ex- 
pected from finite-size effects observed for the KS 
fermion was found for N$ > 16 at (3 = 2.6. It is 
a future problem to see whether and how values 
of Bp vanish toward the infinite volume limit. 

4.2. B K 

The bare value of Br is interpolated as a func- 
tion of m/a with the formula given by chiral per- 
turbation theory, 

Bk = B (1 — 3c(m/<z) log(m/a) -I- 6(m/a)) (5) 

as shown in Fig. |[ The physical value of Bk is 
obtained at half the strange quark mass m s a/2 
(solid circle in Fig. ||) which is estimated from 
the experimental value of mif/m p . 

We plot the renormalized value of 




Figure 4. Scaling behavior of Bk- Previous re- 
sults with the KS action |jj are also shown. 

£>a-(NDR; A* = 2GeV) as a function of the spa- 
tial size in Fig. ||. Filled circles and diamonds are 
data at /3 = 2.6 and 2.9 keeping the same fifth 
dimensional size N$ = 16. The spatial size de- 
pendence is mild especially for N a a > 2.6 fm at 
(3 = 2.6. This result is consistent with that of 
a previous finite spatial size study with the stag- 
gered quark action Q . We conclude that the size 
of about 2.4 fm used in our main runs is sufficient 
to avoid spatial size effects for Bk- 

In Fig. U an open square shows the result for 
N 5 = 32 on a 24 3 x 40 four-dimensional lattice at 
P = 2.6. Since this data is consistent with that 
at A^ 5 = 16 within the statistical error, the fifth 
dimensional size of N$ = 16 is sufficient for Bk- 

Our final results from the main runs are shown 
in Fig. H We observe a good scaling behavior, 
and making a constant extrapolation Bk{o) = 
B K , we find B K (fJ. = 2GeV) = 0.5746(61) where 
the error is statistical only. The open symbols and 
the associated lines represent results from a pre- 
vious calculation with the staggered quark action 
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Figure 5. Renormalized strange quark mass as a 
function of spatial size. 

Q , where gauge invariant and non- invariant four- 
quark operators are used. Compared to the result 
in the continuum limit Bk — 0.628(42) from this 
work, our present result is smaller but consistent 
within the error. 

5. Strange quark mass 

We show results for the strange quark mass in 
the MS scheme at [i = 2 GeV in Fig. ||, where mx 
is used as input. We observe that the fifth dimen- 
sional size N$ = 16 and spatial size N a a = 2.4 fm 
is sufficient for a reliable determination of strange 
quark mass, giving m 8 rj 97 MeV. A recent result 
m s = 110(2) (22) MeV @ using the domain wall 
fcrmions but with the plaquette gauge action at 
(3 = 6.0 is also consistent with ours, if the system- 
atic error from the conversion to the MS scheme 
is taken into account. 

Our results exhibit a good scaling behavior as 
shown in Fig. |(| Making a constant fit yields 
m!p(2GeV) = 97.4(4.4) MeV. Compared with 
results from the conventional 4-dimensional for- 
malisms (the plaquette and Wilson quark action 
fill an d an RG-improved gauge action and the 
clover quark action jll|]), our result is 3-4ct (13- 
18 MeV) smaller. Reduction of systematic errors 
from one-loop renormalization factors and finite- 
size effects would be needed to see if there is an 
actual disagreement. 
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